Complete bond— operator theory of the two— chain spin ladder 
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The discovery of the almost ideal, two-chain spin-ladder material (C5Hi2N)2CuBr4 has once 
again focused attention on this most fundamental problem in low-dimensional quantum magnetism. 
Within the bond-operator framework, three qualitative advances are introduced which extend the 
theory to all finite temperatures and magnetic fields in the gapped regime. This systematic descrip- 
tion permits quantitative and parameter-free experimental comparisons, which are presented for 
the specific heat, and predictions for thermal renormalization of the triplet magnon excitations. 
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I. INTRODUCTION 

Quantum magnetic systems have in recent years of- 
fered a wealth of opportunities for the study of a broad 
range of novel physical phenomena. Systematic improve- 
ments in the growth of pure, single-crystalline samples 
of low-dimensional spin systems, and in the measure- 
ment techniques applied to probe their microscopic prop- 
erties, have resulted in a remarkable convergence of the- 
ory and experiment to explain a number of fundamen- 
tal observations. At the foundation of many of these 
breakthroughs has been the two-chain "spin-ladder" ge- 
ometry of stacked, interacting dimer units, shown in 
Fig. 1. First synthesized in cuprate systems,— both un- 
doped (spin) and doped ladders have been the subject 
of extensive theoretical investigation in the context of 
gapped and gapless quantum magnetism;^ resonating 
valence-bond states?^ generalized frustration effects^i^ 
impurity doping^ (charge) Luttinger-liquid phasesji^ 
and superconductivityi^iiS 

Focusing on dimer-based quantum spin systems, 
examples of physical phenomena realized rather 
well in experiment include magnetic quantum phase 
transitionsiiii^ and quantum critical regimes^ un- 
conventional, critical longitudinal magnonSfi^ the 
Bose-Einstein condensation of magnetic excitationS j^^i^^ 
and the demonstration of their statistical properties J^iii 
However, none of this rich physics is intrinsically a 
consequence of one-dimensionality: while experimental 
efforts to create ladder-like geometries with a range of 
coupling scales and ratios have been rather successful, 
until recently none of the candidate materials has been 
very accurately one-dimensional (ID) in nature. This 
situation has now been changed by the organometallic 
compound (C5Hi2N)2CuBr4, a system of two-chain spin 
ladders whose net interladdcr coupling is approximately 
3% of the characteristic ladder energy scale, which pro- 
pels this geometry back to the forefront of exploration 
into the unique properties of ID magnets. 

(C5Hi2N)2CuBr 4^^'^^ has the additional, major advan- 
tages of forming large, readily deuterated, single crystals 



FIG. 1: (color online) Schematic representation of the spin 
ladder. 5* = 1/2 spins at each lattice site are taken in the 
bond-operator representation to be correlated predominantly 
with the spin on the same ladder rung, forming a singlet state 
(shaded ellipse), except when the rung contains a propagating 
triplet state (depicted as separate spins). J is the exchange 
interaction on a dimer bond, J' the interaction on the chain 
bonds, and / and r are indices for the two chains. The field 
direction sets the triplet quantization axis. 



and of its energy scales falling within the range of labo- 
ratory magnetic fields. The magnetic field is an essential 
parameter in both experimental and theoretical studies. 
The observation, made at a rather early stage, that a spin 
ladder in a field should be driven into a quantum critical, 
spin Luttinger-liquid regime;^ was followed by a number 
of numerical^! and theoretical studie o^^'^'^ of the ground 
state and low-energy excitations. However, the exper- 
imental situation now mandates a physical description 
which captures all of the key quantum fluctuation effects 
in a spin ladder at finite temperatures and applied fields. 

The impressive arsenal of experimental techniques 
available to modern condensed matter physics has in a 
very short time revealed from the (C5Hi2N)2CuBr4 sys- 
tem alone not only most of the phenomena listed in the 
second paragraph but also the field-induced spin Lut- 
tinger liquid. The latter is a direct consequence of the 
impressively ID nature of the material. Thermodynamic 
measurements, specifically of thermal expansion^ spe- 
cific heatr^ and the magnetocaloric effect^i^ have eluci- 
dated the phase diagram: the quantum disordered phase 
gives way to the Luttinger-liquid regime at He = 6.99 
T, full saturation is achieved at Hs = 14.4 T, and 
Luttinger-liquid behavior can be followed up to temper- 
atures of order 1.5Kati7~llT. Nuclear magnetic 
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resonance (NMR)21 has been used to measure critical 
and Luttinger-liquid exponents, the latter varying with 
the applied field. The field-induced 3D ordered phase, 
which is a consequence of the weak interladder coupling 
and sets in below 100 mK, has been found by NMR|22 and 
probed in detail by neutron difFractionj^ Fractionaliza- 
tion of the characteristic ladder magnon excitations into 
field-controlled spinon continua, which constitute unam- 
biguous proof for a Luttinger-liquid description over a 
wide range of energy scales in the lowest triplet branch, 
has been demonstrated by inelastic neutron scattering 
(INS) 3^ 

The spin ladder is the quintessential gapped ID quan- 
tum magnet: with its short correlation length, it is 
perfectly suited for numerical approaches, and many 
of its properties have been computed with great accu- 
racy by exact diagonalization (ED)^ and density-matrix 
renormalization-group (DMRG)^ calculations of the 
ground and excited states, augmented by quantum Monte 
Carlo (QMC}i^ and transfer-matrix rcnormalization- 
group (TMRG)^ studies of thermodynamics. It has also 
been used in the development of systematic, high-order 
techniques such as series expansions^ and Continuous 
Unitary Transformations (CUTs)^ which give quantita- 
tively accurate information concerning high-energy exci- 
tations, including bound states and spectral weights. Si- 
multaneously, however, the spin ladder is rather poorly 
described by most analytical approaches, as much of the 
sophisticated machinery available for ID systems is ap- 
plicable to gapless quantum magnets. Indeed, even the 
numerical approaches listed are all incomplete in some 
regard, most notably concerning finite-temperature dy- 
namics. 

The purpose of this manuscript is to develop one 
straightforward analytical description, based on the 
bond-operator representation of the dimer spin states, to 
give a complete account of the static and dynamic proper- 
ties of the two-chain spin ladder in the spin-gap regime. 
The development requires three qualitative additions to 
the basic bond-operator formalism, specifically finite in- 
terdimer correlation terms, their evolution with field, and 
the incorporation of effective statistics for magnon exci- 
tations. Given the exchange interactions of the ladder, 
the description has no free parameters. 

The focus here will be on the low-field, quantum dis- 
ordered phase of the ladder, where the spin gap is a con- 
sequence of quantum fiuctuations. The high-field, sat- 
urated regime of the ladder is also gapped, and is also 
separated from the Luttinger-liquid phase by a crossover 
(which in higher dimensions is a true quantum phase 
transition), but here the quantum fluctuation effects are 
precisely zero and the gap is a consequence only of the 
field. The zero-temperature properties of this regime 
were described in Ref. [s^. Thermal fiuctuation effects 
do, however, mandate the same treatment of thermo- 
dynamics and finite-temperature dynamics as that pre- 
sented here for the quantum disordered phase. 

The structure of this paper is as follows. Section 



II presents the extended bond-operator formalism. In 
Sec. Ill this is applied at zero temperature to review the 
basic properties and fundamental physics of the spin lad- 
der as a function of the ratio of chain to rung interaction 
parameters. In Sec. IV, the effective statistical proper- 
tics of the magnon excitations are introduced, and em- 
ployed to solve the mean-field equations of the ladder 
at finite temperatures and magnetic fields. Illustrative 
results for the isotropic ladder are used to highlight the 
physics captured by the complete bond-operator descrip- 
tion. In Sec. V, quantitative calculations are performed 
for the parameters of (C5lIi2N)2CuBr4, and particular 
attention is paid to a comparison with experimental data 
for the specific heat. A summary is presented in Sec. VI. 



II. BOND-OPERATOR MEAN-FIELD THEORY 

The Hamiltonian for a Heiscnberg spin ladder in an 
external magnetic field (Fig. 1) is 

H = Y,JSlSl+ ^ J'Sr-S^, + Y.h,-{Sl + Sl) (1) 

i i^m—l.r i 

with h — gfi^H. Here one considers a minimal, purely 
ID model for SU(2) spins with no additional couplings 
(such as frustrating, next-neighbor, three- or four-spin 
interactions) and no anisotropics (which in real materials 
could be of exchange, Dyzaloshinskii-Moriya, single-ion, 
g-tensor, or other origin). The Hamiltonian is trans- 
formed using the identity 

= zL'^^S^ta + ^qs) — i^afs^tj^tj , (2) 

where Sj and ti a{ct = x, y, z) are bond operators for the 
singlet and triplet states of each ladder rungi^i^ These 
operators have bosonic statistics, required to reproduce 
the spin algebra of S\f ^ but they must also obey the local 
hard-core constraint, 

a 

on each ladder rung. This is the physical statement that 
each rung may only be in a singlet state or one of the 
three triplets (equivalent to the four possible states of two 
spin-1/2 entities): these operators arc hard-core bosons. 
In the presence of a magnetic field, and the absence of 
spin-space anisotropy, it is convenient to transform to 
the basis a = 0, — <^ 

In the bond-operator representation, the minimal 
Hamiltonian of Eq. ([T]) takes the form H = HQ + H2+H4, 
where^"— 

Ho - E-J(|4s^ + iO«-)-M»(45.+4,aka-l) 
i 

+ J2Htt+u,+ -tl_U-), (4) 
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with summation over the repeated index a, 

-^2 = +*U4+l,aS»S^+l +H.C. 



with a 



and 



Hi — ^■/'y^[4.o^i+i,o(4.+^-'+i,+ + 



(5) 



(6) 



-t.,ot.+i,o(<+4+i.- + 4,-4+1.+) + H.c] 

The second term in Hq enforces the constraint ([3]) using 
the Lagrange muhiphers /i^. At zero magnetic field, the 
term quadratic in the singlet operators is negative, ensur- 
ing a singlet condensation and justifying the replacement 
Si ~ {si) on each dimer. The ground state of the system 
is then described as a condensate of singlets with a spin 
gap to all triplet excitations. For a uniform ladder there 
is no site variation, and hence one may take (si) = s and 
fii = fj, (at which point it is clear that the constraint is 
enforced only globally, rather than locally). In the ab- 
sence of terms violating the inversion symmetry of the 
rung bonds, there is no three-triplet contribution. The 
four-triplet interaction term is written explicitly to 
facilitate the mean-field decomposition to follow. 

In the gapped phase of the spin ladder, the triplet op- 
erators obey {U^a) = 0, but there is no such condition on 
two-dimer correlation functions. To include the finite- 
field case, one may generalize immediately by introduc- 
ing the three quantities Pa = (t|^ti+i^Q,) and the pair 



i,0''i+l 



(4,+''I+i,-) — (4,-4+1,+) ^ 



-i,o)' Q+- 

Q |_ = Q; in the 0, — basis, all five expectation val- 
ues are real. First introduced in Ref. [HI for the zero- 
field case, the quantities P and Q encode respectively the 
quantum fluctuations for triplet hopping along the ladder 
and for spontaneous creation of triplet pairs on neighbor- 
ing ladder rung bonds. In terms of the spins at each site, 
P and Q describe the nearest-neighbor spin correlation 
functions along the ladder direction, which will be the 
subject of Sec. IIIB. Further, at finite temperature, a 
finite magnetic field will permit a non-zero magnetiza- 
tion 77ii = (t| — {tl_ti,^), which will be uniform 
{nii = m) and whose inclusion constitutes the third and 
final qualitative extension of the basic framework. 

When the Hamiltonian is decoupled using the finite ex- 
pectation values of the preceding paragraph, transformed 
to reciprocal space, and diagonalized in the conventional 
manner;^"— the net result is the quadratic triplet Hamil- 
tonian 



— eo 



J2 ^^."TI. 

k.a 



7fc,a (5 + "(^/c,q)) , (7) 



where the sums for each branch a are normalized to the 
system size and the zeroth-order contribution is 

eo = (-|j-Ai)s2 + |;x-|j-ij'm2 (8) 

P_)-2g2-4QQo] 



-y[Pl + Pl+2Po{P 



per dimer. The operator "fk,a is the bosonic quasiparticle 
diagonalizing the Hamiltonian matrix and the eigenener- 
gics arc 



^fc,+ 



ij' cosA:(P+ -P_ 



J'm 



OJk 



\J' cosk{P+~P-) - J'm 



h, (9) 
(10) 
h, (11) 



with 



Ak 

Afc,o 



Al-Al, 



(12) 



J' cos k (s2 + Po + i(P+ + P_)), 



J' cosfc (s^ - Q - Qo), 

jJ- I-L+ J' COS fc (s2 + P+ + P_), 

J' COS k{s^ - 2Q). 



The quantity n{ujk.a), which enters at finite tempera- 
tures, is the thermal occupation beyond the zero-point 
contribution. In a fully rigorous treatment it is the Bose 
function. Here, however, it will be replaced by an ef- 
fective thermal occupation function which captures cor- 
rectly the hard-core nature of the magnons at high tem- 
peratures. A detailed discussion is deferred to Sec. IV. 

The mean-field equations obtained by taking the 
derivatives of -ffmf © with respect to and s are 

l-s^^y\ ^n{u:k.+) + ^n{iOkfl) + —n[^k,-), (13) 

^ Wfe Ukfl LUk 

k 



J+fi ~ ^ J'cos/c 

k 

. Afe,o- 



Afc-Afc 

n(u}k.- 

^k 



(14) 



Aa:,0 - . Afc — A/. _ 

-I "(w/c.o) H niLOk-) 

UJk UJk 

and to the six finite expectation values are 



Po 

P_ 
Qo 
Q 

m 



cos/e n{uJk,o), 

k ^^'^ 



cos k 



E 



cos k 



^+A\n{uk,+) + (—-l]n{oJk,-) 

^k J \^k 

Afe 



2 

cos k 



l^k 



l)n(t^fc,+) + (^+l)n(tjfe,_) 



cos k 



k 



n[LOkfi), 

l-^k.O 

n[ujk.+ ) H n{ujk-) 

UJk UJk 



^n{ujk,+) - n(wfe,-). 



(15) 
,(16) 
,(17) 

(18) 

(19) 
(20) 



+ 



in which h{uj) = ^ +n{uj). Given the overall energy scale 
J, the self-consistent solution of these equations at each 
value of the field h and the temperature T, which appears 
only in the thermal occupation function, depends solely 
on the parameter ratio J' / J. 
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FIG. 2: (color online) Properties of the two-chain ladder as 
a function of coupling ratio J'/J ai h — — T. (a) Chem- 
ical potential fi. (b) Singlet condensation fraction s = (s). 
(c) Spin gap A: bond-operator data from Eqs. (fT5)) - pU)) 
{BO{c), red, solid line); bond-operator data without P,Q,m 
{BO{2), black, dot-dashed line); ED data for 2x12 system 
{ED 2x12, green circles); ED data extrapolated to an infinite 
ladder {ED, violet diamonds); data from Continuous Unitary 
Transformations {CUT, blue, dashed line), (d) Triplet cor- 
relations P and Q. (e) P and Q for J'/J = 1, shown as 
functions of field at T = 0. (f) Spin correlation ratio B/A 
(see text), comparing bond-operator and CUT results. 



III. ZERO-TEMPERATURE LADDER PHYSICS 

A. Bond— operator description 

An analysis of the physical content of Eqs. ([T^ - (pn|) 
begins by considering the qualitative description of the 
spin ladder at zero temperature and zero field. The pa- 
rameter governing the extent of quantum fluctuations 
away from the limit of decoupled dimers, in which the 
bond-operator approach is exact, but trivial, is the ratio 
J'/J between the exchange interactions on the chain and 
rung bonds. The essential picture is a condensate of rung 
singlets, which are separated by an energy gap of J from 
the (local) rung triplet states. These form a condensate 
in the true sense of the word, its coherence mediated by 
the hopping processes of virtual triplet excitations (the 
quantum fluctuations) along the ladder<^ The effective 
chemical potential fi (Sec. II) is the parameter fixing the 
band center of the triplet magnon modes, and the singlet 



condensate fraction s appears also in the effective band 
width. 

Figure 2(a) shows the evolution of /i as a function of 
J'/J; because the dominant contributions to the integrals 
in the mean-field equations (fT3| -([20 |) are given by those 
k values close to tt, increasingly negative values of /i are 
required to maintain a suitable gap A J' is 

raised (at h = 0, the three modes a are degenerate). The 
singlet condensation fraction s [Fig. 2(b)] decreases from 
precisely unity in the decoupled limit to values of order 
80% beyond the isotropic point (J'/J = 1), suggesting 
a rather high level of singlet condensation even far from 
the strongly dimerized limit. The gap of the spin ladder 
is shown in Fig. 2(c), where the bond-operator technique 
in its simplest, "second-order" forni^ is compared with 
its complete form, and also with data obtained by other 
techniques. ED data for a 2x12 ladder was taken from 
Ref. |39[ and that for the extrapolation to an infinite 
ladder from Ref. [2^. The CUT technique, while per- 
turbative, affords a systematic means of proceeding to 
very high order and can be seen to provide extremely 
accurate results over a wide range of J'/ J for a gapped 
system (i.e. one with a short correlation length) such as 
the two-chain ladder. 

Within the bond-operator description, large quantita- 
tive errors are inevitable at J' > J, where the framework 
loses its validity: the gap begins to rise again as a conse- 
quence of the band center increasing while s falls with J' 
[Fig. 2(b)]. Although the gap prediction is accurate only 
at small J'/J, the second-order approximation delivers 
serendipitously good agreement with the exact result pre- 
cisely at 

while Afj3(l) = 0.501J. 
expectation values P and Q was addressed in Ref. j35l ]. 
a missing factor of 3 for the triplet mode degeneracy 
caused errors in the results for all quantities. In com- 
bination with the above numerical coincidence, this led 
to rather little attention being paid to the P and Q terms 
in the literature. Their qualitative importance will be- 
come clear in Sees. IV and V. Quantitatively, in fact the 
complete bond operator result at the isotropic point is 



the isotropic point, where AgQ 
^ed(I) 



^,'^(1) = 0.503J 
Although the inclusion of the 



A 



3q(1) = 0.640J, a value which is close to half-way 
between the second- and third-order approximations to 
the gap obtained in a strong-dimerization expansion)^ 
Thus one may conclude that the results obtained from 
Eqs. ([T3)) -(|20 |) are essentially optimal within the bond- 
operator framework, and that this is by its nature limited 
to be quantitatively accurate only up to J'/ J ~ 0.5. The 
qualitative physics of the spin ladder may, however, be 
captured for considerably higher values of the coupling 
ratio. 

At finite fields but zero temperature, the three magnon 
bands are split by the Zeeman effect, the uniform linear 
decrease of the entire + branch and increase of the — 
branch with a slope determined by the 5 = 1 nature of 
the magnons and by the T = g-factor. The zero-field 
gap is an important quantity because it determines the 
critical field required to render the lowest mode gaplcss. 
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and thus to drive the system into the Luttinger-Uquid 
regime at the lowest temperatures. 



B. Spin correlations 

Turning to the interdimer correlations P and Q at zero 
field and temperature, the evolution of the two quanti- 
ties as a function of J'/ J is shown in Fig. 2(d). Their 
negative sign is a consequence of the relative phase of the 
triplets. While the triplet hopping term P obviously re- 
mains rather small for all ratios, the pair-creation term 
Q constitutes a very significant quantum fluctuation con- 
tribution even at moderate values of J' / J . For perspec- 
tive, in the isotropic ladder [J' / J = 1) one finds that 
each of the three components [Q appears twice ([72]) ] has 
a value of order 0.23; given that s — 0.83 for J'/J — 1 
[Fig. 2(b)], it is largely the Q terms which exhaust the 
relevant sum rules. The dominance of the triplet pair- 
creation term found here is fully consistent with studies 
of the t-J ladder 

Figure 2(e) illustrates the initially non-intuitive re- 
sult that the zero-temperature triplet correlation func- 
tions do not vary at all with the applied magnetic field: 
P_i_ = Pq = P_ and Q = Qo at all fields, reflecting the 
unbroken SU(2) spin symmetry of the electronic wave 
function in the gapped phase. However, it becomes clear 
by inspection of the mean-field equations (fT3 |) - (f20l) that 
in fact a finite temperature is required before the applied 
field can affect either the P and Q terms or the magne- 
tization TO. This behavior is the topic of Sec. IV. 

There is a direct relationship between the quantities P 
and Q and the interdimer spin correlations (5^-S''_|_j) = 
= B: from Eq. © it follows that, up to 
fourth-order corrections, B = |s^(-P + Q)- Further, 
the on-dimer spin correlation (S^-S^) = A is straight- 
forwardly A = which clearly has the required 
limit as J'/J and s — > 1, and hence the ratio 
B/A = -2(P + Q). This ratio has been measured^iii 
by neutron-scattering studies of (C5Hi2N)2CuBr4, and 
the results are consistent with the ratio J'/J ~ 1/3.9 de- 
duced for this material (Sec. V). Bond-operator results 
for B/A are compared in Fig. 2(f) with the CUT analy- 
sis of the spin ladder, also presented in Ref. [iOl, over a 
range of coupling ratios J'/J. As for the gap, the bond- 
operator technique delivers quantitative accuracy up to 
J'/J ~ 0.5, while deviations appear and increase beyond 
this. 

Thus the intrinsic nature of the quantum mechanical 
wave function of the spin ladder, by which is meant the 
extent of spin correlations across the different bonds, may 
be deduced directly from the triplet correlations. Their 
evolution due to quantum fluctuations (this Sec.) and 
to thermal fluctuations (Sec. IV), as well as the com- 
bined effects of both types of fluctuation in a magnetic 
field, may be followed in the bond-operator approach. 
This type of analysis may be applied to a number of 



gapped, low-dimensional quantum spin systems, includ- 
ing the dimcrized = 1/2 chain Cu(N03)2-2.5D20i^ and 
the 2D 5 = 1/2 system PHCC^ for which systematic 
measurements of the quantities Jij{Si ■ Sj) have already 
been performed. 

IV. FINITE TEMPERATURES 

Turning now to the situation at finite temperature, 
there is naturally no explicit expression for a thermal oc- 
cupation function of hard-core bosons. This problem was 
addressed for the ladder system in Ref. [3 by enforcing 
the local constraint ([3]) at the global level to obtain an ef- 
fective partition function for the hard-core boson system. 
The effective statistics deduced from this expression, 

"K,a,/5) = — ^F^— 7^, (21) 

1 +z;a^"(/3) 

where (3 is the inverse temperature and Za[l3) = 
X]fe e~''"'='° is the partition function of mode a, were 
shown in Ref. [13 to give a spectacularly successful de- 
scription of the thermal band-narrowing and gap en- 
hancement measured in the 3D couplcd-dimer system 
TlCuCls. However, it should be emphasized that the 
global reweighting of the partition function formulated 
for a spin ladder in Ref. [l6j is entirely dimension- 
independent, and thus should function equally well for 
all hard-core boson systems. Indeed the same ansatz has 
also been used to obtain good descriptions of thermody- 
namic quantities in the 2D dimcrized spin-chain system 
(VO)2P207>^ and of spectral weights for the excitations 
of ultra-cold bosonic atoms on a ID optical lattice4^ 

It is important to recall the meaning of Eqs. ([13])- 
(|20p at finite temperatures. The spectrum of excitations 
of any Hamiltonian is a temperature-independent quan- 
tity, and the temperature is responsible only for driving 
changes in the weights with which given modes are occu- 
pied. These spectral weights are represented by the dy- 
namic structure factor S{k, uj, T), which can be measured 
by INS. The bond-operator approach delivers an effective 
description of the one-magnon sector of the dynamical 
structure factor, in the form of a peak position for each 
of the (field-split) magnon branches. Changes in these 
peak positions as a function of temperature are in fact a 
statement of the leading thermal effects on the spectral- 
weight distribution: the key result of Ref. [T^] is that 
thermal shifts in the peak positions by factors of up to 
3 are not accompanied by thermal broadening of the ex- 
citations over the same energy width. That the magnon 
modes remain coherent shows that their spectral-weight 
shifts are a true many-body effect governed by the overall 
band structure and by their hard-core statistics. These 
features are captured over the full range of temperatures 
(up to T J) in the bond-operator framework with 
Eq. (EI]). 

Hcuristically, a band-narrowing effect arises because, 
while the temperature does not alter the magnetic ex- 
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FIG. 3: (color online) (a) Magnon bands Uaik), shown for 
the isotropic ladder at T/J = 0.1 and in a field h/J = 0.5. 
(b) Thermal evolution of the band maximum, center, and 
minimum for a spin ladder with J'/J = 1 and h = Q. 



change interactions, thermal occupation of some dimer 
sites by triplets acts to block the motion of a propagat- 
ing test triplet due to the local hard-core constraint (O . 
This is reflected in the fact that the effective occupation 
function for mode a (|2ip involves explicitly terms due 
to all three magnon branches, thus encoding a form of 
inter-branch exclusion. In the self-consistent solution of 
the mean- field equations, an increase in n{uJk.a) with in- 
creasing temperature causes a drop in all of the quantities 
s, P, and Q. The physical interpretation of this effect 
is a suppression both of on-dimer singlet condensation 
and in the coherence of interdimer fluctuation processes, 
which results in a narrowing of the magnon bands and 
consequent increase in the gap. While such an effect is 
also obtained with bosonic statistics, only the effective 
occupation function of Eq. (|2T|) takes into account cor- 
rectly the number of available hard-core boson states of 
the system. 

In this section, the qualitative physics of the two-chain 
spin ladder at finite fields and temperatures is illustrated 
using the isotropic coupling ratio, J'/ J = 1. This value is 
not chosen for any special reason: it is already clear from 
the previous section that the bond-operator description 
is not quantitatively accurate beyond J'/J = 0.5, and it 
is not the case that some of the qualitative effects become 
stronger as the coupling ratio increases (indeed, as shown 
in Sec. V, the converse is true for certain quantities). The 
choice is simply a good representative value for a ladder 
away from the strongly dimerized limit, and serves as 
an instructive counterpoint to the results presented in 
Sec. V for J'/J = 1/3.9 (a value which is effectively in 
this limit). 



A. Magnon dispersion 

Figure 3(a) serves first as a reminder of the magnon 
band structure: there are three branches, which all have a 
minimum at the antiferromagnetic wave vector k = n and 
which are similar in shape and splitting. However, for a 
ladder beyond the strongly dimerized limit (J' <C J), and 
at finite fields and temperatures, these branches are not 
symmetrical (in k) about a band center, they do not have 
exactly the same line shape, and their separation is not 
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FIG. 4: (color online) (a) Triplet correlations P (blue, dashed) 
and Q (red, solid) as functions of temperature for J'/J — 1 
and h — 0. P is shown multiplied by a factor of 4 to illustrate 
its relative thermal renormalization. (b) Triplet correlations 
P+, Po, and P- shown as functions of field at T/J = 0.1. (c) 
P+, Po, and P_ and (d) Q and Qo as functions of temperature 
for J'/J = 1 and h/J = 0.5. 



exactly the bare Zeeman splitting defined in Sec. IIIA; 
all of these facts can be read from Eqs. ([9|)- p2|) . 

Considering first the case with no magnetic field. 
Fig. 3(b) illustrates the thermal band-narrowing effect 
by showing for the isotropic ladder the T-dependence of 
the band maximum (fc = 0), band minimum (fc = tt), 
and band center defined as the mode energy at fc = 7r/2. 
The suppression of hopping begins when T becomes sim- 
ilar to the gap, and the effect saturates when T exceeds 
J. The form of the thermal renormalization, and of 
the consequent growth of the gap with temperature, is 
exactly analogous to the results obtained for TlCuCla 
in Ref. [l3|. As noted above, uJk,a obtained from the 
bond-operator formalism is the peak position at finite 
temperture of the one-magnon spectral weight of the 
spin ladder, a quantity not yet available from any nu- 
merical approaches. However, some recent attempts to 
compute effective magnon spectra for quantum magnets 
at finite temperatures, specifically ED calculations for 
short, dimerized (gapped) S" = 1/2 chains^ and real- 
time DMRG calculations for gapless S" = 1/2 Heisenberg 
chains;^ suggest that further progress can be expected 
in this direction. 

The evolution of the zero-field triplet correlations P 
and Q corresponding to this band-narrowing is shown 
in Fig. 4(a). The two quantities have contrasting be- 
havior: while finite temperatures act only to suppress 
pair-creation processes (i.e. rendering these less coher- 
ent), the hopping correlations show an initial rise, or 
thermal enhancement of such processes, before thermal 
decoherence becomes the dominant effect. Here the terms 
"coherence" and "decoherence" are used in the sense of 
Sec. HI, the former to indicate the intrinsic correlations 



FIG. 5: (color online) (a) Band-narrowing effects for the 
isotropic ladder at h/J = 0.5. (b) Gap-curvature effects for 
the band minimum of the three triplet branches, illustrated 
for the isotropic ladder with T/J = 0.1. 



driven by quantum fluctuations (which depend on J'/ J, 
Fig. 2) and the latter to indicate the suppression of these 
correlations by thermal fluctuations. 

At finite fields and temperatures, the full complexity of 
the system becomes visible, with effects appearing which 
arc not present at finite h ov T alone. Figures 4(c) and 
(d) show the thermal evolution of the triplet correlation 
functions at finite h: P-^ is strongly and Q weakly en- 
hanced at small fields, before thermal decoherence takes 
over; P_ is generally anticorrelated to P+, but has no 
mechanism for initial suppression corresponding to the 
enhancement of P+ ; Pq shows a non-monotonic evolution 
as a consequence of the competing effects of coupling to 
P-i_, the field scale, and the two opposing forms of ther- 
mally driven behavior [Fig. 4(a)]; Qo shows the same 
weak anticorrelation to Q, followed by thermal decoher- 
ence [Fig. 4(d)]. The correlation functions vary some- 
what predictably with h when T is fixed, as shown in 
Fig. 4(b), although it should be noted that the flat line is 
P_, while Pq is the function which drops more strongly 
as P+ is enhanced. This result can be understood by 
considering the corresponding temperature in Fig. 4(c). 

It is important to note at this point that it is not par- 
ticularly meaningful to discuss temperatures T in excess 
of, or even on the order of, the rung coupling scale J. 
At these temperatures, the concept of the rung singlet 
and triplets as well-defined states of the two rung spins 
breaks down, and it is more appropriate to consider these 
spins separately. Even at temperatures considerably be- 
low J, magnon excitations become broad functions in the 
space of energy and wave vector: the complex question 
of this intrinsic line width in a gapped quantum mag- 
net has been considered only at the scmiclassical level^ 
and more microscopically in some limiting cases^ but 
is not addressed within the bond-operator framework at 
the current level of approximation. Thus the data up 
to T = 2J shown in Figs. 3 and 4 are to be taken only 
as an example of the physics encoded in the mean-field 
equations, and all of the figures to follow will be cut at 
T = J. For those concerning the critical properties of the 
ladder near the field-driven crossover to the Luttingcr- 
liquid phase, further caveats will be necessary. 

The thermal rcnormalization of the triplet band width 
at finite field is shown in Fig. 5(a). The effect of the 



FIG. 6: Magnetization m as a function of T for four different 
magnetic fields (a) and as a function of h for three different 
temperatures (b). 

field is to reduce further the symmetry of the band edges, 
causing an initial upward trend at the upper edge as ther- 
mal fiuctuations allow the polarizing effect of the field to 
operate. Analogous to the effects visible in Figs. 3 and 
4, this behavior is reversed at higher temperatures, as T 
competes with h, resulting in an overall downward trend 
of the band center. The consequences for the mode gaps 
as a function of field at finite T are shown in Fig. 5(b): 
here one may regard the 0-mode curvature as an indica- 
tion of generic thermal effects (increased net energy), on 
top of which the additional contributions contained in the 
P and TO terms in Eqs. (j9]) and (fTTj) act against this cur- 
vature for the upper-mode gap, giving a rather straight 
field-dependence, but an enhanced curvature away from 
gap closure of the lower mode. This reflects the basic 
physical phenomenon that the lower mode absorbs most 
of the thermally excited magnons at flnite T and h, and 
is thus the most strongly rcnormalized; the necessary de- 
viation from Zeeman splitting is contained in the P and 
TO terms. 



B. Thermodynamic properties 

It is instructive to consider the magnetization to to 
understand its evolution with fleld and temperature, and 
hence how it affects the magnon dispersion branches. For 
a flxcd field as a function of temperature [Fig. 6(a)], to is 
non-monotonic, first rising as a consequence of the ther- 
mally excited states which become available to react to 
the field, but becoming maximal as T approaches J. The 
low-temperature rise is clearly of activated form, the gap 
diminishing as the field is raised. Only close to the criti- 
cal field which closes the gap {h/J = 0.64 in the isotropic 
ladder) does the exponential form turn over to a power 
law, which in ID is expected to be a square root as a 
function of temperature at very low T. The value of m 
at the maximum is clearly an almost linear function of the 
field which is polarizing the available triplet states. The 
high-P behavior, found also in extensive TMRG studies 
of m{h,T)^ sets in when T competes with the inter- 
action energy scales of the system, suppressing not only 
triplet propagation and the polarized moment, but also 
dimer singlet formation itself. 

Figure 6(b) shows m as a function of h for different 
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temperatures. In general, both the number of states be- 
coming available and the field which is able to polarize 
them give contributions to m. At low T the response is 
activated and the dominant contribution is due to the 
field. At high T, for which T/ J = 0.5 is already sufla- 
cient, m becomes an entirely linear function of h: here the 
temperature ensures a large number of available states, 
and these are polarized by whatever field is applied. This 
physics, also illustrated below in Figs. 9(e) and (f), con- 
stitutes a regime in which the gap renormalization [of 
the type shown in Fig. 5(b)] is also completely linear. 
This results in a "pseudo-Zeeman" behavior of the mode 
minima in an applied field, where the gaps evolve lin- 
early but with a slope reduced from the conventional, 
zero-temperature value by the P+ — P- and m terms 
in Eqs. (O and (HJ). It is possible that this behav- 
ior has been misinterpreted in some previous analyses 
of low-dimensional quantum magnets, for example as a 
T-dependcnt g-factor. 

At low temperatures one finds that m ~ \{P- ^ P+)i 
and thus that the two contributions to non-Zceman be- 
havior of the gap are essentially equal. This value of m 
can be regarded as the quantum-fluctuation contribution 
to the magnetization. As T is increased, m increases 
compared to |P-f — P-|, until at high temperatures it 
is exponentially greater: while m is enhanced by ther- 
mal population effects, there is no such contribution to 
|P+ — P-l, which indeed is suppressed at high tempera- 
tures [Fig. 4(c)]. 

In dimensions of 2 and higher, it is of great interest 
to investigate the temperature-dependence of the criti- 
cal field he- the field which closes the excitation gap of 
the quantum spin system, driving a quantum phase tran- 
sition. The quantity hc{T) — hc{0) encapsulates the crit- 
ical properties of the theory, and has been analyzed in 
detail for the higher-dimensional cases, where the tran- 
sition is in the Bose-Einstein universality class In this 
case, the transition is from a (singlet-based) quantum 
disordered state to an ordered magnetic one, the order 
parameter being the transverse (staggered) moment. The 
low-temperature behavior of the phase-transition line 
should have the scaling form hc{T) - hc{0) = AT'^^^, 
where A is a constant, d is the spatial dimensionality of 
the system, and the denominator 2 in the exponent en- 
ters from the quadratic dispersion of the free hard-core 
bosons {lo (X k^) at the transition. 

In a ID system, there is no ordered phase or or- 
der parameter, and indeed no transition: hc{T) marks 
a crossover to the Luttingcr-liquid regimci^i^ Under 
these circumstances, the definition of a quantity hc{T) 
may not be unique. The crossover has nevertheless been 
found to have a very clear signature both in experimen- 
tal measurements of the thermodynamic properties of 
(C5Hi2N)2CuBr4 and in numerical studies of the magne- 
tization m{h, T). In measurements of the magnetocaloric 
eSect^ its location was deduced from the minima and 
maxima in m{T) through the condition dm/dT\h = 0, 
while in TMRG calculations^ it was taken from the 




FIG. 7: (a) Critical magnetic field required to close the 
effective spin gap of ladders with coupling ratios J' / J = 1 
and 0.25 at finite temperatures, (b) h^°{T) - h^°{0) on 
logarithmic axes to illustrate square-root scaling (solid lines). 

maximum in dm/dh. The quantity h^^{T) investigated 
in the bond-operator description is that value of the 
field where the maximum of the finite-temperature, one- 
magnon spectral function is moved to zero frequency, 
and therefore does not necessarily have a direct relation 
to the crossover field (which is dictated by the lowest 
level in the excitation spectrum): it is rather a state- 
ment about the reaction of the spectral-weight function 
to thermal flucutations, which may be very different in 
ID from the higher-dimensional cases. For practical pur- 
poses, h^'~'{T) marks the limit of applicability of the 
bond-operator formalism at the level of Sec. II, which 
cannot be used for the gapless phase without extensive 
modification. At T = 0, the quantity he = h^'~'{Q) is 
the same crossover field as that obtained from all other 
definitions. 

The critical field h^'~'{T) is shown in Fig. 7(a) for the 
isotropic ladder, and by way of comparison also for a 
ladder with coupling ratio J'/J = 1/4. At low temper- 
atures, both curves have a very rapid increase with T: 
the logarithmic fit in Fig. 7(b) reveals this singular be- 
havior to have precisely the square-root form obtained 
from the analysis above with d = 1. Thus the effect 
of thermal fiuctuations on the low-temperature spectral 
function is indeed dramatic. Beyond the critical regime, 
/i^'-'(r) becomes fiatter at higher temperatures, meaning 
that the maximum of the excitation spectrum becomes 
rather temperature-insensitive. A slight decrease is visi- 
ble for J'/J = 1/4 as the system approaches the "break- 
down" regime T ^ J, where the theory contains dimin- 
ishing energy scales [visible in the magnon bands (Figs. 3 
and 5) and in the magnetization (Fig. 6)]. Returning to 
the critical regime, while it may be taken to exist over a 
temperature range of order 0.2 J in the isotropic ladder, it 
is clearly correspondingly narrower for smaller coupling 
ratios (T/J - 0.05 at J'/J =1/4). This narrow ID crit- 
ical regime is precisely that part cut ofi' by 3D coupling 
even in very weakly interacting, quasi-lD ladder systems 
such as (C5Hi2N)2CuBr4. 

A further thermodynamic quantity which can be com- 
puted in this framework is the magnetic specific heat. 
Obtained from the second temperature derivative of the 
effective hard-core boson free energy, this can be ex- 
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pressed in the form 



(22) 



Consideration of the magnetic specific heat is deferred to 
Sec. V, where it is discussed in the context of a direct 
experimental comparison. 



V. (C5Hi2N)2CuBr4 

Having understood in full the physics contained in 
the complete mean-field equations, one may turn to the 
nearly ideal spin-ladder material (C5Hi2N)2CuBr4. Fits 
to the magnon bands, of the type shown in Fig. 3(a), 
were performed in Ref. [2^ and will not be repeated here. 
In the bond-operator framework they give J = 13.07 
K, J' = 3.35 K, and hence a coupling ratio within the 
spin ladders of J'/ J = l/3.9<^ In this strongly dimerized 
regime (sometimes referred to as the "strong-coupling" 
regime), the bond-operator description is extremely ac- 
curate (Sec. III). 

The real material has a number of complications which 
may affect the accuracy of a detailed comparison between 
this theoretical description and experimental measure- 
ment on levels of up to 5%. For the present purposes, 
the most important is the anisotropy in the g-factor) ^^i^° 
which enters the comparisons between thermodynamic 
measurements (most performed with H \\ a) and INS 
dispersion relations [H \\ b). The value chosen here to 
reproduce He = 6.99 T, g ~ 2.12, is close to the aver- 
age in the (a6)-plane of the system. Among the other 
materials questions, one is the weak interladder coupling 
mentioned in Sec. I, which is approximately 3% of 
another is a detectable magnetostriction effect (changes 
of J and J' with applied field) of order 1%, which be- 
comes relevant when comparing the low- and high-field 
gapped regimes?^ a third is a spin anisotropy, or a break- 
ing of SU(2) spin symmetry of as-yet undetermined ori- 
gin, which a recent ESR investigation^ has estimated as 
high as 5% of J. While a spin anisotropy can be included 
in the current framework when enough information be- 
comes available, the high-field regime is not the focus of 
the current manuscript and the 3D nature of the material 
is not of prime importance while the system has a robust 
spin gap, meaning at all fields below where h enters the 
critical regime close to he- 

In this section, results are presented only for ID lad- 
ders with coupling ratio J'/ J = 1/3.9. Fig. 8 shows the 
triplet correlations P and Q which assist a qualitative un- 
derstanding of the ladders in (C5Hi2N)2CuBr4. Figure 
9 shows the magnetization and one-magnon band struc- 
ture of these ladders at finite fields and temperatures, 
giving direct experimental predictions. Figure 10 con- 
siders the thermodynamics of these ladders through the 
magnetic specific heat, including a direct experimental 
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FIG. 8: (color online) Triplet correlation functions for the 
ladders of (C5Hi2N)2CuBr4 at different temperatures and ap- 
plied magnetic fields, (a) Pa for a = +,0, — as functions of 
temperature at a field of 1 T. (b) Qo and (5+- = Q as func- 
tions of temperature at fields of 1 T and 6 T. (c) Pa at a field 
of 6 T. (d) Pa as a function of field for temperatures of 1 K 
and 5 K. 



comparison and quantitative reproduction of published 
data for this quantity. 

Beginning with the triplet correlations, the quanti- 
ties and Qa. are shown in Fig. 8 for ladders with 
J' I J = 1/3.9. It is clear in Fig. 8(a) that the general 
enhancement of hopping correlations by a finite temper- 
ature is large in strongly dimerized ladders, and that the 
scale for thermal decoherence to dominate the correlation 
physics of (C5Hi2N)2CuBr4 is T ~ 10 K. The same quan- 
tity shown not at low field but close to the critical point 
[Fig. 8(c)] emphasizes the very strong effect which the 
field can have on the enhancement of in this regime 
of J' I J [note the different scales in panels (a) and (c)]. 
The triplet pair-creation correlation function Qq, how- 
ever, depends strongly on neither field nor temperature, 
being sensitive only to the thermal decoherence effect 
[Fig. 8(b)]. As a function of the magnetic field [Fig. 8(d)], 
splitting of the hopping correlations in the three branches 
requires significant field strengths at low temperatures, 
where thermal occupations are activated, but the effect 
is nevertheless very strong on approaching the critical re- 
gion. By contrast, at higher temperatures where states 
in all branches are available, all applied fields are capable 
of splitting the branches but the effect is rather modest. 

Figure 9 turns to some of the obscrvablcs of the ladder 
system. The magnetization as a function of tempera- 
ture over a range of fixed magnetic fields is shown in 
Fig. 9(a): initial activated behavior dictated by the spin 
gap leads to a peak around T ~ J/2 and then to the 
slow decrease indicative of the loss of coherence in the 
magnetic sector. These features are qualitatively similar 
to those of the magnetic specific heat, shown in Fig. 10, 
but differ in magnitude because m increases as the sys- 
tem approaches the critical field [Fig. 6(a)]. Figure 9(b) 
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FIG. 9: (color online) Calculations for (C5Hi2N)2CuBr4 lad- 
ders at finite temperature and magnetic field: as functions of 
temperature, (a) magnetization m for fields of 1 T (black, 
solid line), 3 T (green, long-dashed line, 5 T (blue, dot- 
dashed line), and 6.99 T (red, dashed line) and (b) magnon 
band width at a field of 2 T; as functions of the applied mag- 
netic field, magnetization m (c,e) and gaps (d,f) at tem- 
peratures of 1 K (c,d) and 10 K (e,f). 



shows the thermal band-narrowing effect of Figs. 3(b) 
and 5(a), illustrated for a magnetic field of 2 T. This 
phenomenon was explained in Sec. IV, and INS measure- 
ments are expected shortly. Data for the magnetization 
m and for the gap of each mode as functions of the 
applied field, at temperatures equivalent to 1 K and 10 
K, are shown in Figs. 9(c)-(f). The magnetization at 1 K 
[Fig. 9(c)] shows the activated behavior visible at lower 
temperatures in Fig. 6(b), and the gap of the low-lying 
mode [Fig. 9(d)] has a significant curvature of the type 
shown in Fig. 5(b). By contrast, the 10 K data are well in 
the high-temperature limit discussed in Sec. IV, where 
the ready availability of thermally excited states leads 
to a linear magnetization [Fig. 9(e)] and to the pseudo- 
Zecman form of the mode gaps as functions of tempera- 
ture [Fig. 9(f)]. It should be stressed, however, that in 
experiment it is somewhat unlikely that clearly defined 
magnon excitations would still be visible at temperatures 
as high as 10 K, due to extensive thermal broadening of 
the response function. If one considers temperatures of 
approximately 5 K, where a definite magnon mode may 
still be detectable, then these data are rather similar to 
the forms shown in Figs. 9(e) and (f), albeit with some 
discernible curvature remaining in m and A_|-. 
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FIG. 10: (color online) (a) Magnetic specific heat of 
(C5Hi2N)2CuBr4 at magnetic fields of T (black, uppermost 
peak), 3 T (green, second peak), 5 T (blue, third peak) and 
7 T (red, lowest peak). Solid lines are given by the complete 
bond-operator theory, while crosses show experimental data 
presented in Ref. [23]. (b) Specific-heat data from bond- 
operator theory, shown on logarithmic axes to illustrate the 
approach to a low-T square-root form as h ^ he- The red, 
solid line marks a slope of 1/2. 



The most directly measurable and perhaps informa- 
tive single thermodynamic quantity is the magnetic 
specific heat. Detailed measurements performed on 
(C5Hi2N)2CuBr4 over the full range of temperatures, 
and of applied magnetic fields in both the gapped (quan- 
tum disordered) and gapless (Luttingcr-liquid) phases, 
are presented in Ref. [25|. The experimental results 
were compared with somewhat numerically demanding 
ED and DMRG studies of ladders with the same cou- 
pling ratio, and also to a bond-operator approximation 
in the gapped regime and a Bethe-Ansatz approxima- 
tion in the gapless one. The same behavior of the mag- 
netic specific heat as a function of field and temperature 
has also been demonstrated in TMRG studies of ladder 
thermodynamics,— but has not been computed explicitly 
for the parameters of (C5Hi2N)2CuBr4. 

The magnetic specific heat given by the complete 
bond-operator approach is shown in Fig. 10(a) for fields 
corresponding to T, 3 T, 5 T, and 7 T, along with ex- 
perimental measurements made at the same fields.— In 
the regime with a robust spin gap, the specific heat shows 
exponential activation at low temperatures, followed by 
a peaking at T ~ J/2. The activation gap decreases as 
the field is raised, resulting in more weight at lower tem- 
peratures and a drop in the peak height. By contrast, 
the critical field which closes the triplet gap in a single 
(C5lIi2N)2CuBr4 ladder is 6.99 T, and hence the sys- 
tem shows gapless behavior at 7 T; note that while the 
bond-operator description is no longer valid at h > he 
for T — 0, the rapid rise of h^'-' with temperature shown 
in Fig. 7 ensures that the gapless nature of the excita- 
tions becomes a numerical problem only at temperatures 
two orders of magnitude below those shown in Fig. 10(a). 
As the field approaches he, the peak height at T ^ J/2 
continues to fall; both this peak and the high-T decrease 
are features sensitive not to the gap but to the center and 
upper edges of the excitation band. 

The low-T form of the specific heat in the gapless 
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regime is a power law, and the power in a spin ladder at 
the critical field is 1/2. Figure 10(b) shows specific-heat 
data obtained as h approaches he from below: clearly a 
square-root dependence is obtained, but only for fields 
extremely close to the critical value (c/. Ref. [2J|). That 
even differences of 0.01 T can have a significant impact 
on the form of thermodynamic quantities indicates at the 
theoretical level the possibility that the quantum critical 
regime may be extremely narrow, and at the experimen- 
tal level the need for exquisite accuracy in the measure- 
ment of critical properties. 

The data computed from Eqs. ([T^ - ipO)) . represented 
as lines in Fig. 10(a), is in perfect quantitative agreement 
with essentially every feature of the and 3 T experi- 
mental data, and shows only the most minor deviations 
at 5 T. There is no separate scale factor for each curve, 
i.e. as already stressed above there are no free parame- 
ters in the theory. Only on the critical line, at the very 
limit of applicability of the bond-operator framework, 
are some small but systematic discrepancies noticeable: 
specifically, the complete bond-operator result gives a 
peak height (around 5 K) which is marginally too low, 
and this missing weight appears to be present around 
temperatures of 2 K. However, the results of Fig. 10(b) 
indicate that extreme caution is required before ascribing 
this to a failure of the bond-operator description, as the 
true proximity of the experimental measurements to the 
critical field is accurate perhaps only to within 0.05 T. 

The bond-operator approximation applied in Ref. (25l | 
used the same effective magnon statistics as here, but 
not the P, Q, and m terms; the same qualitative form of 
error was found when h ^ he 'as in Fig. 10(a), but it was 
very much larger (of order 15% in the peak height). Fig- 
ure 10(a) makes clear that the self-consistent inclusion 
of these terms results in a very significant correction of 
this error. Thus the complete bond-operator framework 
is highly accurate even at the very limit of its regime of 
applicability and even in one dimension. It was noted 
in Ref. [25} that some degree of error may be introduced 
by the treatment of the hard-core constraint ([3]): while 
the total number of states is correct, global application 
of the constraint may result in some of these states not 
being penalized sufficiently for the fact that they violate 
the constraint at the local level. Extra weight would then 
be expected not in the peak but at lower energies, per- 
haps resulting in the apparent persistence of the low-T 
critical (power-law) regime to artificially high tempera- 
tures. However, the results in Fig. 10 show that these 
errors are at the percent level even for h = he, and that 
the complete bond-operator theory is capable of a fully 
quantititively accurate description of all features of the 
gapped, strongly dimerized spin ladder. 

A more explicit summary statement on the accuracy of 
the complete bond-operator treatment for the two-chain 
spin ladder, drawn from all of the results above, is the fol- 
lowing: it is extremely high, meaning quantitative agree- 
ment within a few percent when compared to numerical 
results (particularly DMRG for the spectrum and TMRG 



for thermodynamics) for the strongly dimerized regime 
J ^ J'. This level of agreement, which is achieved with 
no adjustable parameters, extends to J' ~ J/2 and ap- 
plies to all physical observables computed here (spin cor- 
relations, magnon spectra, thermodynamic quantities). 
Quantitative deviations from percent levels of accuracy 
become noticeable as J' /J is further increased towards 
the isotropic case, a result not at all surprising in view of 
the expected regime of applicability of the bond-operator 
method, and the type of quantum fluctuations which it 
is designed to capture. In the isotropic ladder, where 
the qualitative physics of the system is nevertheless very 
well described within the bond-operator approach (hence 
its use for illustrative purposes in Sec. IV), the level of 
quantitative error as judged from the interdimer spin cor- 
relations and the zero-field spin gap is of order 20%. 



VI. SUMMARY 

A complete version of the bond-operator theory for the 
spin ladder has been developed, which captures all of the 
relevant physics of the quantum disordered system, at all 
temperatures and all magnetic fields up to the closure of 
the spin gap. [At this point the magnon excitations frac- 
tionalize into spinons and a Luttinger-liquid description 
is appropriate.] At the heart of the finite-temperature 
framework is a global ansatz for the effective exclusion 
statistics of the magnon excitations, which are hard-core 
bosons. Among the phenomena described with quan- 
titative accuracy and with no free parameters are the 
spin correlations, the thermal evolution of the magnon 
bands, and their thermodynamic properties. Extensive 
comparison becomes possible with existing data for the 
spin-ladder material (C5Hi2N)2CuBr4, and predictions 
are made for forthcoming experiments which will investi- 
gate the thermal renormalization of the magnon excita- 
tions. The same treatment is also applicable to calculate 
finite-temperature properties in the quantum disordered 
phases of other couplcd-dimer systems of any dimension- 
ality. 



Acknowledgments 

The authors are grateful to T. Giamarchi, M. Sigrist, 
and G. S. Uhrig for helpful discussions, to E. Gull for 
computing assistance, and to K. P. Schmidt and G. S. 
Uhrig for provision of the CUT data shown in Figs. 2(c) 
and (f). The Pauli Center for Theoretical Studies at ETH 
Zurich and the Department of Physics of the University 
of Fribourg extended their generous hospitality to BN 
during the completion of the manuscript. This work was 
supported by the National Science Foundation of China 
under Grant No. 10874244, by Chinese National Basic 
Research Project No. 2007CB925001, and by the Royal 
Society. 



12 



^ Z. Hiroi, M. Azuma, M. Takano, and Y. Bando, J. Solid 

State Chem. 95, 230 (1991). 
^ see E. Dagotto and T. M. Rice, Science 271, 618 (1996), 

and references therein. 
^ S. R. White, R. M. Noack, and D. J. Scalapino, Phys. Rev. 

Lett. 73, 886 (1994). 

* N. Okazaki, K. Okamoto, and T. Sakai, J. Phys. Soc. Jpn. 
69, 2419 (2000). 

^ V. Gritsev, B. Normand, and D. Baeriswyl, Phys. Rev. B 

69, 094431 (2004). 
® M. Sigrist and A. Furusaki, J. Phys. Soc. Jpn. 65, 2385 

(1996). 

C. A. Hayward and D. Poilblanc, Pliys. Rev. B 53, 11721 
(1996). 

* H. H. Lin, L. Balents, and M. P. A. Fisher, Pliys. Rev. B 
56, 6569 (1997). 

^ E. Dagotto, J. Riera, and D. Scalapino, Phys. Rev. B 45, 
5744 (1992). 

^° M. Sigrist, T. M. Rice, and F. C. Zhang, Pliys. Rev. B 49, 
12058 (1994). 

see T. Giamarchi, Ch. Riiegg, and O. Tchernyshyov, Na- 
ture Physics 4, 198 (2008), and references therein. 
Ch. Riiegg, A. Furrer, D. Sheptyakov, Th. Strassle, K. W. 
Kramer, H.-U. Giidel, and L. Melesi, Phys. Rev. Lett. 93, 
257201 (2004). 

" Ch. Riiegg, B. Normand, M. Matsumoto, A. Furrer, D. Mc- 
Morrow, K. Kramer, H.-U. Giidel, S. Gvasaliya, H. Mutka, 
and M. Boehm, Phys. Rev. Lett. 100, 205701 (2008). 
T. Nikuni, M. Oshikawa, A. Oosawa, and H. Tanaka, Phys. 
Rev. Lett. 84, 5868 (2000). 

Ch. Riiegg, N. Cavadini, A. Furrer, H.-U. Giidel, K. 
Kramer, H. Mutka, A. Wildes, K. Habicht, and P. Vorder- 
wisch. Nature 423, 62 (2003). 

M. Troyer, H. Tsunetsugu, and D. Wiirtz, Phys. Rev. B 
50, 13515 (1994). 

Ch. Riiegg, B. Normand, M. Matsumoto, C. Niedermayer, 

A. Furrer, K. Kramer, H.-U. Giidel, Ph. Bourges, Y. Sidis, 
and H. Mutka, Phys. Rev. Lett. 95, 267201 (2005). 

B. R. Patyal, B. L. Scott, and R. D. Willett, Phys. Rev. B 
41, 1657 (1990). 

1^ B. C. Watson, V. N. Kotov, M. W. Meisel, D. W. Hall, 

G. E. Granroth, W. T. Montfrooij, S. E. Nagler, D. A. 

Jensen, R. Backov, M. A. Petruska, G. E. Fanucci, and D. 

R. Talham, Phys. Rev. Lett. 86, 5168 (2001). 
^° S. Sachdev, T. Senthil, and R. Shankar, Phys. Rev. B 50, 

258 (1994). 

C. A. Hayward, D. Poilblanc, and L. P. Levy, Phys. Rev. 
B 54, R12649 (1996). 

^2 R. Chitra and T. Giamarchi, Phys. Rev. B 55, 5816 (1997). 
F. Mila, Eur. Phys. J. B 6, 201 (1998). 
T. Lorenz, O. Heyer, M. Garst, F. Anfuso, A. Rosch, Ch. 
Riiegg, and K. Kramer, Phys. Rev. Lett. 100, 067208 
(2008). F. Anfuso, M. Garst, A. Rosch, O. Heyer, T. 
Lorenz, Ch. Riiegg, and K. Kramer, Phys. Rev. B 77, 
235113 (2008). 

Ch. Riiegg, K. Kiefer, B. Thielemann, D. F. McMorrow, 
V. Zapf, B. Normand, M. B. Zvonarev, P. Bouillot, C. 
KoUath, T. Giamarchi, S. Capponi, D. Poilblanc, D. Biner, 
and K. W. Kramer, Phys. Rev. Lett. 101, 247202 (2008). 
B. Thielemann, Ch. Riiegg, K. Kiefer, H. M. R0nnow, B. 



Normand, P. Bouillot, C. KoUath, E. Orignac, R. Citro, 
T. Giamarchi, A. M. Lauchh, D. Biner, K. W. Kramer, F. 
Wolff-Fabris, V. Zapf, M. Jaime, J. Stahn, N. B. Chris- 
tensen, B. Grenier, D. F. McMorrow, and J. Mesot, Phys. 
Rev. B 79, 020408(R) (2009). 

M. Klanjsek, H. Mayaffre, C. Berthier, M. Horvatic, B. 
Chiari, O. Piovesana, P. Bouillot, C. KoUath, E. Orignac, 
R. Citro, and T. Giamarchi, Phys. Rev. Lett. 101, 137207 
(2008). 

B. Thielemann, Ch. Riiegg, H. M. R0nnow, J.-S. Caux, A. 
M. LauchU, B. Normand, D. Biner, K. W. Kramer, H.-U. 
Giidel, J. Stahn, S. N. Gvasaliya, K. Habicht, K. Kiefer, 
M. Boehm, D. F. McMorrow, and J. Mesot, Phys. Rev. 
Lett. 102, 107204 (2009). 

^® T. Barnes, E. Dagotto, J. Riera, and E. S. Swanson, Phys. 

Rev. B 47, 3196 (1993). 
^° X. Wang and L. Yu, Phys. Rev. Lett. 84, 5399 (2000). 

S. Trebst, H. Monien, C. J. Hamer, Z. Weihong, and R. R. 

P. Singh, Phys. Rev. Lett. 85, 4373 (2000). 

C. Knetter, K. P. Schmidt, M. Griininger, and G. S. Uhrig, 
Phys. Rev. Lett. 87, 167204 (2001). 

B. Normand, Acta Phys. Pol. B 31, 3005 (2000). 
S. Sachdev and R. Bhatt, Phys. Rev. B 41, 9323 (1990). 
S. Gopalan, T. M. Rice, and M. Sigrist, Phys. Rev. B 49, 
8901 (1994). 

M. Matsumoto, B. Normand, T. M. Rice, and M. Sigrist, 
Phys. Rev. Lett. 89, 077203 (2002); Phys. Rev. B 69, 
054423 (2004). 

^'^ B. Normand and T. M. Rice, Phys. Rev. B 54, 7180 (1996). 
B. Normand and T. M. Rice, Phys. Rev. B 56, 8760 (1997). 
M. Reigrotzki, H. Tsunetsugu, and T. M. Rice, J. Phys. 
Cond. Matter 6, 9236 (1994). 

A. T. Savici, G. E. Granroth, C. L. Broholm, D. M. Pa- 
jerowski, C. M. Brown, D. R. Talham, M. W. Meisel, K. 
P. Schmidt, G. S. Uhrig, and S. E. Nagler, Phys. Rev. B 
80, 094411 (2009). 

B. Thielemann, Ph. D. thesis, ETH Zurich (2009). 

G. Xu, C. Broholm, D. H. Reich, and M. A. Adams, Phys. 
Rev. Lett. 84, 4465 (2000). 

M. B. Stone, L Zaliznyak, D. H. Reich, and C. Broholm, 
Phys. Rev. B 64, 144405 (2001). 

G. S. Uhrig and B. Normand, Phys. Rev. B 58, R14705 

(1998); Phys. Rev. B 63, 134418 (2001). 
'^^ A. Reischl, K. P. Schmidt, and G. S. Uhrig, Phys. Rev. A, 

72, 063609 (2005). 
■"^ H.-J. Mikeska and C. Luckmann, Phys. Rev. B 73, 184426 

(2006). 

''■^ T. Barthel, U. SchoUwock, and S. R. White, Phys. Rev. B 
79, 245101 (2009). 

K. Damle and S. Sachdev, Phys. Rev. B 57, 8307 (1998). 
F. H. L. Essler and R. M. Konik, Phys. Rev. B 78, 
100403 (2008); A. J. A. James, F. H. L. Essler, and 
R. M. Konik, Phys. Rev. B 78, 094411 (2008); W. D. 
Goetze, U. Karahasanovic, and F. H. L. Essler, unpub- 
lished (arXiv:1005.0492). 

E. Cizmar, M. Ozerov, J. Wosnitza, B. Thielemann, K. 
W. Kramer, Ch. Riiegg, O. Piovesana, M. Klanjsek, M. 
Horvatic, C. Berthier, and S. A. Zvyagin, unpublished 
(arXiv:1005:1474). 



